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Abstract 

We have diagonalized the transfer matrix of the U q [osp(2\2m)] vertex model by means 
of the algebraic Bethe ansatz method for a variety of grading possibilities. This allowed 
us to investigate the thermodynamic limit as well as the finite size properties of the 
corresponding spin chain in the massless regime. The leading behaviour of the finite size 
corrections to the spectrum is conjectured for arbitrary m. For m = 1 we find a critical 
line with central charge c = — 1 whose exponents vary continuously with the g-deformation 
parameter. For m > 2 the finite size term related to the conformal anomaly depends on 
the anisotropy which indicates a multicritical behaviour typical of loop models. 
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1 Introduction 



Two-dimensional vertex models of statistical mechanics are nowadays considered classical 
paradigms of the theory of exactly solvable models [I]. Their statistical weights can be di- 
rectly related to the elements of a i?-matrix satisfying the Yang-Baxter equation invariant 
relative to the fundamental representations of U q [Q] quantum symmetries [2]. 

The thermodynamic limit properties of most vertex models derived from ordinary Lie al- 
gebras, such as the free-energy and the nature of the excitations, have been well examined 
over the past decades in the literature, see for instance [3], SJ [5j [6] and references therein. It 
is believed, for instance, that the massless regimes of these vertex models are described by 
the critical properties of Wess-Zumino-Witten field theories on the group Q [7J. We remark, 
however, that at least one counter example to such common belief appears to occur in the 
U q [sp(2m)] vertex models [Sj- 

By way of contrast, similar physical properties of the U g [Q] vertex models when Q is a 
superalgebra have not yet been examined in details. The majority of the results concerning 
the possible universality classes of critical behaviour governing the massless phases in these 
systems have been concentrated on the sl(n\m) symmetry [9j QUI EH]. Similar information for 
other superalgebras such as osp(r\2m) has so far been restricted to the rational limit q — ► 1 
[T2| [T3l [H] . It is not yet clear, however, if the determined classes of universality are robust 
against (/-deformations such as the cases of ungraded algebras. 

In this paper we hope to start to bridge this gap by investigating the leading finite size 
corrections governing the eigenspectrum of the U q [osp(2\2m)} vertex models. These finite size 
properties have a direct relationship with the critical operator content of massless phases |15j . 
In order to do that we have diagonalized the respective row-to-row transfer matrix by means of 
the algebraic Bethe ansatz approach. We have considered explicitly all grading choices that are 
compatible with the underlying U(l) symmetries of the i?-matrix. This step will complement 
our previous efforts concerning the Bethe ansatz solution of the U q [osp(n\2m)] vertex models 
|16j . We recall that the exact solution for n = 2 was not presented before [16J due to technical 
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problems with the special grading considered in that work. Here we are able to circumvent 
such technicalities. 

This paper is organized as follows. We start next section by describing the statistical 
weights of the U q [osp(2\2m)] vertex models. In section 3 we discuss the diagonalization of 
the corresponding row-to-row transfer matrix, within the algebraic Bethe ansatz method, for 
a variety of grading possibilities. In section 4 we use such grading freedom to choose the 
appropriate one in order to deal with the thermodynamic limit in the simplest possible manner. 
In section 5 we study the finite size properties of the U g [osp(2\2m)] vertex models by both 
analytical and numerical approaches. This provides us the basis to conjecture, in the massless 
regime, the behaviour of the leading finite size corrections to the spectrum for general m. For 
m — 1 these results indicate that the central charge of the underlying conformal field theory 
is c = — 1. For m > 2, however, we find that the finite size term associated to the conformal 
central charge depends on the anisotropy coupling q. In Appendices A and B we describe the 
technical details entering the Bethe ansatz solution of a particular grading. 



2 The U q [osp(2\2m)] vertex model 

The i?-matrix of the U q [osp(2\2m)] vertex model is defined on the tensor product of Z 2 graded 
spaces having two species of bosons and 2m species of fermions. The Grassmann parity p a is 
used to distinguish the bosonic p a = and the fermionic p a = 1 degrees of freedom. 

To establish the statistical interpretation of this system it is important to know the structure 
of the i?-matrix i? 12 (A) in appropriate coordinates such as the Weyl basis. This task is in 
general rather involved for superalgebras but recently some progresses towards this direction 
have been made [TTJ [18] . The Boltzmann weights of such systems can be conveniently written 
in terms of the standard relation [19J, 

fli 2 (A)=P 12 i? 12 (A), (1) 

N 

where Pi 2 is the graded permutator given by P = (— l) PaP0 e a p ® ep a and e a p denotes 

a,/3=l 
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N x N matrices having only one non-null element with value 1 at row a and column {3. The 
operator ^12 (A) satisfies the following form of the Yang-Baxter equation, 



£12 (A - »)Rx{X)RM = j R 23 (m)^i2(A) j R 23 (A - /i). 



(2) 



which is insensitive to grading. 

It turns out that the corresponding ^-matrix of the U q [osp(2\2m)] vertex model, in terms 
of the Weyl basis, can be written as, 



N 



N 



a = l 



a,/3=l 



N No 

+c (0) (A) e Qa ® e^ + c (0) (A) £ e aQ ® 



a,/3=l 



a, 0=1 
a>j3,a^f3> 



No 
a,/8=l 



(3) 



For later convenience we have introduced the label (0) = (2|2m). It emphasizes that we 
are considering a Z 2 graded space with two bosonic and 2m fermionic degrees of freedom 
and Nq — 2 + 2m denotes the dimension of such space. Each index a has its conjugated 
a' = N + 1 — a and the Boltzmann weights a^'(X), tf® (A), c^(A) and c^(A) are given by 

a (°)(A) = ( e 2 >-C( ))(e 2A ( 1 -^ 0) )-g 2 e 2A ^ 0) ) 6(°)(A) = g(e 2A -l)(e 2A -C (0) ) 
C (°>(A) = (l- g 2 )(e 2A -C (0) ) C (°)(A)=e 2 V°)(A), 



(4) 



while d^g(A) has the form 



f g(e 2A - l)(e 2A - C (0) ) + e 2A (g 2 - 1)(C (0) - 1) 

2A 



<S(A) 



;e- - 1 

,2 



(g 2 - 1) 

(g 2 - l)e 2A 



e 2A_ c (0) )( _ 1)P L 0) g 2 p (°) +e2 A (g2 _ 1) 

C(°)( e 2A -l)^^-^(e 2A -C (0) ) 

e /3 



,2A 



-i)-g^-^(e 2A -C (0) ^ 



e /3 



a = /3 = /3 > 
a = [3^/3' 

a < (3 
a > (3 



(5) 



We stress that the formulas ([MS]) are valid only for grading choices whose respective parities 



p^ satisfy the reflexion condition = p ( "/. These grading possibilities are consonant with 
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,(°) 



the underlying U(l) symmetries of the system that usually play an essential role in Bethe 
ansatz solutions. Furthermore, the parameter = q~ 2m and the variables e a and t a are 
related to the parities by 



„«>) 



1 < a < 



iVn 



p (0) 
-1) 2 



(6) 



+ 1 < a < N 



a 



a 



1 _ D (o) 

2 la 



2 E $ 



(o) 



a</3<- 



(0) 



1 < a < f 1 



(7) 



^ + l<a<N 



"^ 0) + 2 E 

^+l</3<a 

We would like to close this section with the following remark. The above explicit expression 
for the U q [osp(2\2m)] .R- matrix was first presented by us for a particular grading choice (16] and 
later on generalized to include other grading possibilities satisfying the condition = p^) |17j . 
In the former reference we claimed also to have exhibited the explicit expression of the ^-matrix 
associated to the twisted U q [osp^ (2n\2m)\ quantum superalgebra. Recently, however, we 
realized that such identification is not correct and the ^-matrix denoted by U q [osp^ 2 \2n\2m)] 
in [12] is in fact the one invariant relative to the U q [spo(2n\2m)] quantum symmetry Q [2D]. This 
means that the results to be obtained in next sections are therefore also valid for the vertex 
model based on the U q [spo(2m\2)] symmetry. We believe that the correct U q [osp^ (2n\2m)] 
i?-matrix were indeed obtained by us in [17] as those associated with the generalizations of 

(2) 

Jimbo's -D„+! Z?-matrix. We hope that this later identification could be confirmed in near 
future by means a detailed analysis of the set of algebraic relations coming from the respective 
Yang-Baxter algebra [21]. We also note that the i?-matrices associated to g-deformations of 
the osp(2\2) symmetry have been previously investigated in 



1 We thank J.R. Links for suggesting us that this may be the case. 
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3 The algebraic Bethe ansatz 



The quantum inverse scattering method provides us a systematic framework to construct and 
solve integrable vertex models by the algebraic Bethe ansatz [21]. It also can be extended to 
systems whose i?-matrices are invariant relative to Lie superalgebras [19]. In this approach we 
start by considering a collection of i?-matrices, R_Aj(X) with j — 1, . . . , L, acting non-trivially 

L 

on the auxiliary space = C^ and on the j-th. node of the quantum space C^ . An 
important ingredient is the monodromy matrix defined by the following ordered product of 
i?-matrices, 

T(°)(A) = i?S(A)4°Li(A) • • • i2g(A). (8) 

The row-to-row transfer matrix of the respective vertex model can then be written as the 
supertrace of the monodromy matrix with respect to the auxiliary space [19], namely 

N ° (0) 

T(°)(A) = Str^o)[r(°)(A)] = £(-l) pL) TW(A). (9) 

a=l 

The next step is to present the solution of the eigenvalue problem, 

T(°)(A)|$) =A (0) (A)|$), (10) 

within an algebraic formulation of the Bethe ansatz. 

In this section we tackle the problem ffTUl) in the case of the i?-matrices ([3]17j) of previous 
section for any of the grading p^ = p^) choices. We remark that such solution for a variety 
of such gradings is in general rather intricate even for the U q [sl(n\m)] vertex model [251 [2S]- 
Here follow the nested Bethe ansatz formalism developed in [27] for isotropic vertex models 
and recently extended to accommodate trigonometric i?-matrices based on g-deformed Lie 
superalgebras [15] • We recall, however, that in the later reference the Bethe ansatz solution 
for the specific case of the U q [osp(2\2m)\ vertex model was not presented and here we will be 
filling this gap. Considering that the main procedure has already been well explained before 
[27J EES] there is no need to repeat it again in details. In what follows we shall restrict ourselves 
only to the essential points concerning the solution of such eigenvalue problem. Fortunately, 
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we find that the presence of the many grading possibilities p^ = p^) can still be accommodate 
in terms of certain recurrence relations envisaged by us in [TB] for a specific grading choice. 
This relation for the eigenvalues of T^(X) turns out to be, 



A^)(A,{AS a \...,A£)}) = (-i)^n(-i)^ 4 q) (a-a^) n(-i)^ 1 { i 



fc(a)( A f +1 )-A) 



+ (-i)« n(-D«« A (A - a<« >) n (-^ £ t Hx ~\Ju 

»=1 i=1 d N a ,N a { A ~ A i ) 



+ Il^ ) (A-A l H) n , a) '"\ (a+1) A^^A^A^,-,^)}), 

i=l »=1 «V Q ,V Q l A _ A i J 

(11) 

while the corresponding Bethe ansatz equations for the rapidities {A^ } are given by, 

i=l ) ijtj d N a ,N a { A j ~ A i > a l \ A i ~~ A j ) 

n a+ i , (a+1) / ■> (a+1) \( a+1 h 

x TT ( i) p i 1 ^ - j j ri2) 

L > ^ + i)( A ^)-Af +1) ) • ( j 

We now describe the way the recurrence relations (I11|I12I) should be interpreted. The label 
(a) in the eigenvalues and Bethe ansatz equations was introduced to characterize the respective 
graded space these results are concerned with. The dimension of such space is twice less than 
the one we started with, N a = Nq — 2a, and the respective number of bosonic and fermionic 
degrees of freedom are determined by the following rule, 

(a) = (tf a -5>£>|5>W). (13) 

(3=1 13=1 

The Grassmann parities pp associated with the graded space (a) are obtained through 
the relation Pg" = Pb+i f° r P = 1j • • • > N a — 2. The Boltzmann weights a[ (A), (A) and 
Na (\) are derived from (J3HZD, considering the graded space characterized by (a) instead 
of the original one labeled (0). Finally = (— l) p i ) q 1 ~ 2p i ' and the consistency with the 
original eigenvalue problem requires us to set Xj — for j — 1, . . . , no and to make the 
identification n = L. 



In order to obtain the eigenvalues and respective Bethe ansatz equations for a given choice 
of parities p$ = p^) we need to iterate the relations (11 111 121) starting from a = 0. We then 
carry on such nested procedure until we reach a final step labeled by (/) and therefore up to 
a — f — 1. In this last step we have to deal with the diagonalization of an inhomogeneous 
transfer matrix of the following type, 



r(/)(A,{A^,...,AW}) 



Str 



Ad) 



R 



>(/) 



?(/) 



(14) 



The solution of the eigenvalue problem for such last step depends much on the choice of the 
parities we started with. We find that for all gradings choices satisfying p^ = p^), except the 
special case p^ = 1 for a — 1, . . . , m, m + 3, . . . , 2m + 2 and p^+i — Pm+2 = 0, the last step 
consists in the diagonalization of a common six-vertex model. In our notation it is identified 
as (/) = (0|2) and the i?-matrix governing such final step is, 

fa[ f \X) 



v o 



i? (/) (A) 



with the following Boltzmann weights 

a[ f \\) = (e 2A - ? - 4 )(e 2 V-l) 



42(A) = ^(g 4 -l)(e 2 V-l) 






\ 




iS (a) 42(A) 







a (A) 42(A) 





■ 





4 /} (A)J 




42(A) = £ 


(e 2A -l)(e 2 V 


42(A) = \ 

Q 


e 2X (q 4 - 


l)(e 2A 



(15) 



(16) 



Considering that the Bethe ansatz solution of the six vertex model has been already well 
examined in the literature, we shall not extend over this problem. In order to present our 

^ (a) 

results in a more suitable form we define Q a (X) = sinh (A — X\ ) and set q = e n . In this 

8=1 

way we have the following expression for the eigenvalues 



A(°)(A) 



;_1)p 



(0) 



(-i)^ 0) 4 0) (A) 



l Qi (A + (-l)^f) 

Qi (A-(-iH 0) z|) 



(0) 

+ (-l)^o 



+ 



1) P "°4U(A) 
^ (0) (A)]'E^(A|{Af}) 



Qi (A+(2m-(-l)fi 0) )i|) 
Qx (A + (2m+(-lH 0) ) ? f) 



2m 



a=l 



(17) 



where the auxiliary functions G Q (A|{A^}) are given by, 



Ga(A|{Af}) 



= < 



(0) (0) 

t (A-^-(-l) p "+i» 7 ) Q a+1 (\-8 a+1 + (-lf<*+in) 



' Qa(A-5a) Q a+ l(A— <5 a+ l) 

Qm(A-<5 m +i7) Qm+i(A-5 m +i-2t7) 

Qm(A-lSm) Qm+l(A — <5 m + l) 



a = 1, . . . , m — 1 

a = m 

a — m + 1, . . . , 2m 



[ G Q _ m Hm 7 - A| - {Af }) 
The rapidities |A^| are constrained to satisfy the following set of Bethe ansatz equations, 



n Q _i sinh A 



n 



i sinh(Af ) -Ar i) + (-l)^° ) z|) 



|| sinh (xf - X^ + z£; q7 



x 



+j sinh (aJ q) - X\ a) - ik a ^j 

sinh (a^ - Xf - (-l)^^l) 



2 

jot 2 



nm sinh fAS m+1) -Al m) +i 7 
II 



5 sinh (aS q+1) - Xf + (-iy«U 9o 

a — 1, ... ,m 

"2+1 sinh (V m+1) - A, (m+1) + 2 ?7 



= A i sinh (Aj m+1) - A, (m) - i 7 ) 



n 



^ sinh (Af +1) - Af +1) - 2* 7 ) 



(18) 



where k n = 



(-1)^ + (-1)^ 



and # 



a = m 



I 1 otherwise 

We remark that in order to obtain the Bethe ansatz equation in the above symmetric form 
we have performed the shifts {A^} — > {X^} + 8 a . The variables 5 a have a strong dependence 
on the parities and are given by, 

„(°) 



^E(-i)^ 



• 7 
*2 



f3=l 
m 

E(-i) 3 



(0) 



a — 1, . . . , m 
a = m + 1 



(19) 



As usual we see that the Bethe ansatz equations as well as the eigenvalues depend strongly 
on choice of the parities . This feature has been captured here in a unified way by the index 
k a . The possible different forms of Bethe ansatz equations concerning the distinct grading 
choices can be better appreciated in terms of Dynkin diagrams. In this representation the 
scattering factors between the rapidities and Xj are recasted in terms of the elements 
e aj 3 of the respective Cartan matrix. In order to be more specific we exhibit in Figure 1 the 
diagram related to the grading 



a = 1, N, 



pL o) = t ' (20) 

l 1 otherwise 



l 



m m+1 



Figure 1: Representation of the Bethe ansatz equations (fT5|) in the grading BF . . . FF . . . B. 
The other grading possibilities in the family considered so far, namely 



f a = 3,3' 3>1 
P^ = \ , (21) 

I 1 otherwise 



are represented in Figure 2. 




(3-1 B m m+1 



Figure 2: Representation of Bethe ansatz equations (JTSJ in the grading F ... FBF . . . FF . . . FBF ...F. 

We now turn our attention to the Bethe ansatz solution for the remaining grading, 

, n , f a = m + 1, m + 2 
PL 0) = • (22) 

I 1 otherwise 

For the grading choice ( 1221) the last step is no longer governed by the six-vertex model. 
In this last stage one has to deal with a 16 x 16 i?-matrix R^'{X) whose graded space is 



(/) = (2 1 2). This problem is in fact a special case of the one associated with the general 
U q [osp(2\2)] i?-matrix built from the admissible one-parameter four-dimensional representa- 
tion [23]. The Bethe ansatz solution of such vertex model in the grading (1221) involves extra 
technicalities such as the presence of auxiliary transfer matrices that cannot be written as 
trace of monodromy operators. Here we avoid overcrowding this section with more technical 
details and we summarized them in Appendices A and B. In order to solve the nested problem 
for the U q [osp(2\2m)] vertex model one needs to use the final results given in Eqs.( 1A.25llA.26l) 
together with the recurrence relations (lll)12p . By performing these steps we find that the 
corresponding eigenvalues in the grading (122]) are, 

, 2m 

A(°)(A) = Ao(A) - [6<°>(A)] E^(A|{Af}) 

a=l 



Ao(A) 



-4 0) 



-4 0) (A) 



(A) 



d {0) 
a N ,N 



L Q+(\-q) g-fA-*j) 



(A) 



Qi(A+(2m+l)i; 
Q 1 (A+(2m-l)ii) 



A+«; 



a+«; 



-^Svivo(A) 



q + (a+^| 7 ) 



(A+i| 7 ) 



m > 1 
m = 1 



(23) 



G a (X\{\f}) 

■ Q a {\+( a +2)q) Q a+1 (A+(q-l)jl) 
Q a (\+a.q) Q a+ i(\+(a+l)il) 
Q,n-i(X+(m+l)q) Q+(\+(m-2)i%) Q_ (A+(m-2)i%) 
Q m -l[\+(m-l)i^ Q+^A+rm^ 
Q+(X+(m-2)q) Q_(A+(m+2)i£ 



)- 1 A+mi^ 



] + ( X+mi'-. 



Q-{\+mq) 



[G Q _ m (-zm 7 -A|-{Af ) }) 



a 

a 

a ■ 
a - 



l,...,m-2 
m — 1 



m 

m+ 1, 



2m 



provided the rapidities {A; } satisfy the following Bethe ansatz equations, 
>-i) 



n «-i sinh f A^ a) 

n 1 



A 



+ q 



i=i 



sinh (A 



(a) 



A 



(a-l) 



1 2 



Hm-2 

n 



sinh (A 



(m-2) _ \(m-l) .7 



sinh (Xj 



m-2) 



-A 



(m-1) _ -7 
4 2 



S sinhfA^ + 
i^j sinh (A^ q) 



A • - «7 



n 

n TO _i 

n 



sinh (a, {q+1) - Af + q) 

sinh (A 
sinh ^A 



(a+l) 
(m-l) 



Af - z| 

(rn— 1) 



- A 



# sinh(Aj m 1} 
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A 



(m-l) 



a 



m-2 



n+ 



n 



sinh (A 



(+) 



i sinh ( A 



(+) 



yim-l) 



+ q) ^ sinh (Aj~ 
sinh (\- 



n 

i=i 



A 



(m— 1) 



+ <i) 



(m-1) 



sinh (Af } - Al m - 1} 



n 

i=l 



sinh (\j 



±) 



A 



(m-1) 



i 1 



n 



sinh f 




-Af> 


-27) 


sinh ( 




-Af» 


+ 27) 



(24) 



We recall that the symmetrical form ( 1241) is obtained after performing the shifts {Aj } 



1 m 



1, and {A^ } — > {A^} — im^. We close this section by 



A±) 



presenting in Figure 3 the diagrammatic representation of the Bethe ansatz equations ([2 

+ 




Figure 3: Representation of Bethe ansatz equations (J24J) in the grading F ... FBBF ...F. 



4 Thermodynamic Limit 

In this section we will study the thermodynamic limit properties of the spin chain associated 
to the U q [osp(2\2m)} vertex model. The corresponding Hamiltonian is formally obtained as 
the logarithmic derivative of the transfer matrix (Q at the regular point A = 0, 

H = -J[T^]- 1 ^-T^(\)\ X=0 , (25) 
aA 

where from now on we have fixed the normalization J = sin(7). 

We start our analysis by studying the spectrum of the operator (1251) for small chains, in 
the anti-ferromagnetic regime < 7 < |, by means of exact diagonalization methods for 
m = 1,2. The next step is to reproduce the lowest energies within the Bethe ansatz solutions 
of previous section in order to find the pattern of the corresponding roots {X^}. This helps 
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us to select, among the possible forms of the Bethe ansatz equations, the one that has the 
less complicated root structure as possible. This study leads us to select the set of Bethe 
ansatz solution associated to the grading F . . . FBBF . . . F, since the low-lying spectrum of 
(1251) is reproduced by using mainly real roots for all nested levels. In this case, the eigenvalues 
i? (m) (L, 7) of the Hamiltonian ([25}, up to an additive constant, are given in terms of the 
variables {A^ } by, 



-2 J sin (7) 



E^(a! 1} ,7) m>\ 



' 1 (26) 



n+ n_ 



i=i i=i 



where e(A, 7 ) = cosh(2A) _ cos(7) ■ 

We now explore the Bethe ansatz equations on the grading F . . . FBBF ... F in order to 
determine analytically the ground state energy and the nature of the low-energy excitations. 
Considering that the low-lying spectrum is described mostly in terms of real roots we take 
directly the logarithmic of the original Bethe ansatz equations ff24l) and as result we find, 

5 ltl L*{\f>l)=2*Q? + £$(A? ) -A« 7)- £ £*(A?-Ai Q) ^) 

Z fc=l a=J±lfc=l 

Z = l,...,m-2 (27) 
S^mxf^^Qf + -^,7)- Y. E*(Af-Ar»,|) 

Z fc=l a=m-2,± fe=l Z 

/ = m - 1 (28) 

S m ,m>!!\ 2) = 2 T Qf + £ *(A<" - Af , 7 ) - *(A«" - Af-", 2) 

Z fc=l fe=l Z 

/ = ± (29) 

where $(A, 7) = 2 arctan [cot (7) tanh (A)]. The numbers Qj define the different branches of 
the logarithm and in general are integers or half-integers. For example, part of the low-lying 
spectrum can be parameterized in terms of an integer sector index r\ and the corresponding 
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sequence of Q)p numbers are 



Q 



(i) 



Q) 



[L-n-i]+j-i 



L 



— r± — 1 



L — T\ and 1 = 1, 
3 = 1, 



,m-l (30) 
L 



- r± 



(31) 



For large L, the number of roots n\ tends toward a continuous distribution in the real axis 
and the following density of roots can be defined 

1 



p«(A«)= lim — , ( . 

L -*°° L (Xj+i — Xj ^ 



(32) 



In the limit L — > oo, the Bethe ansatz equations (I27H29P turn into coupled linear integral 
equations for the densities p^(X^). These integral equations can be solved by standard Fourier 

/+°° 
e iU * x p">(X)dX, and the final results are 
-oo 

lCOsh -fc#l 1 = 1 m-1 

(33) 



2tt cosh HHii 



47r cosh [2nii] 



1, 
± 



From Eqs. ( l26|33l) we can compute the ground state energy per site e^\^) in the infinite 
volume limit. By writing Eq. ff26l) . in terms of its Fourier transform, we find the expression 



3 (m) 



(7) 



-2J 



sinh 




cosh 


(m - 


-l)f 


sinh 


TTLJ 

2 


cosh 







do; 



(34) 



Let us now turn our attention to the behaviour of the low-lying excitations in the ther- 
modynamic limit. As usual to many integrable models, the energy e^ l '(x) and the momenta 
(x) of the l-th. excitation measured from the ground state are related by 



2np {l \x 



p {l \x) 



(y)dy 



(35) 



By using Eq. (1331) we conclude that the low-momenta dispersion relation is linear for all 
the excitations, e^(p) = Vm^p"'. The respective sound velocity v m {^) is found to be 



V m {l) = — ■ 



(36) 



We have now the basic ingredients to investigate the finite size effects in the spectrum of 
the ?7g[(2|2m)] spin chain for < 7 < |. 
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5 Finite size properties 



The basic behaviour of the leading finite size corrections to the spectrum of gapless systems 
are expected to follow that of conformally invariant theories in a strip of width L [15]. For 
periodic boundary conditions, the ground state energy E (L) behaves, for large L, as 

m. = 6oo - ™ + { L-\ (37) 

where c is the central charge and v is the sound velocity. 

The structure of the higher energy states E a (L) are also determined by the conformal 
dimensions X a of the respective primary operators, namely 

E a (L) E (L) 2nvX a 



L L L 2 



0{L~ A ). (38) 



In what follows we begin our study of the finite size effects by considering first the simplest 
case m = 1. 



5.1 The U q [osp{2\2)} model 



The finite size corrections for the C/g[osp(2|2)] spin chain can be studied with rather little effort 
at the particular point 7 = §• In this case the Bethe ansatz equations for the roots {A^- } 
become similar to that of lattice free-fermion models, 

L 



sinh (xf } + %\ 



iLk 



(±) 



;_!)«, 



sinh (Af ) - %\ 
while the spectrum are parameterized by 

n-\. n_ 

EW(L,-) = -2 J ^ cos (A:j +) ) -2J^cos (fcH 



(39) 



(40) 



i=i j=i 

Therefore, for the value 7 = f , one can exhibit exact expressions for the low-lying energies, 
in a given sector r±, by summing over selected free-momenta of type = j^nf where nf are 
integers. The computation depends, however, whether the fermionic index r + + r_ is an odd 
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or an even number. When r + + r_ is an odd number we find that 

-2J 



7Y . 



whose asymptotic expansion for large L becomes 



cos 




+ cos (^) 




sin 







(41) 



,7I\ 



27T 



,7T, 



1 

6 



r 2 _|_ r 2 
4 



+ L 



(42) 



Similar analysis can be performed when r + + r_ is an even number. The difference is that 
now the free-momenta are shifted by a fixed amount -|. For instance, the expression for the 



lowest energy in the sector r + 
£«(L,f) 



is given by 



-4Jtan (— % 



3 (i) 



,7I\ 



,7]\ 



27T 

Z^ ll 2 



1 
6 



O (L 



-2 



(43) 



Direct comparison between Eqs. fl42|43l) reveals us that the form of the finite size effects 
has a clear dependence on the fermionic index r + + r_. We also note that the ground state 
for finite L lies in the odd sectors r + = ±1 and r_ = or r + = and r_ = ±1 and therefore 
it is four-fold degenerated. This preliminary analysis will be of utility to help us to make a 
prediction for the finite size corrections in the whole regime < 7 < |. 

To make further progress for arbitrary values of the coupling 7 we use the so-called density 
root method [28l [29] . This approach is able to give us the main expected behaviour of the 
leading finite size corrections when both the ground state and the low-lying excitations are 
described in terms of real roots or those carrying a fixed imaginary part such as if. This 
is exactly the situation we found for the U q [osp(2\2)] model. This conclusion is achieved by 
comparing the spectrum generated by numerical solutions of the Bethe ansatz equations (1241) 
with that from direct diagonalization of the Hamiltonian (12"51) up to L = 16. By applying the 
root density approach to the [7 g [osp(2|2)] model one finds that its prediction for the finite size 
behaviour of the eigenenergies is 



£«(L,7) 



2ty 



^()) + ^i(7) 



__ + x s +' s ~ 



(7) 



+ o (l- 2 ) , 



(44) 
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where the scaling dimensions X*+'*z(l) can be written as 



7\( r + + r -) 2 7 (r + — r_) 2 1 (s + + s_) 2 7r (s. 



\2 



+ ■ 7T 4 7T 4 (1 — ^) 4 7 4 

The integers r± are related to the number of roots m± — \ — t± while the numbers s± are 
directly related to the presence of holes in the distribution. The latter indices are rather 
sensitive to boundary conditions and therefore they need extra care. In fact, by comparing 
Eqs. (I44II45I) at the point 7 = | with Eqs. ( 14"2H4"5|) one sees that for r + + r_ odd the number 
s± are expected to start from zero. By way of contrast for r + = r_ = the lowest values for 
s± is in fact half- integer |. From our numerical analysis we also conclude that the standard 
root density assumptions concerning the values for s± are valid only for anti-periodic boundary 
conditions in the case r + + r_ is an even number. This means that in such sectors, integers 
values for s± are expected only when a twist e ±47r multiplies the Bethe ansatz equations (1241) for 
both {A^} variables. It turns out that the effect of a twist e l<fl in the root density method is 
to shift the numbers s± by a factor Considering these observations one concludes that, for 
periodic boundary conditions, the numbers s± should indeed begin at values ±| when r + + r_ 
is even. These arguments strongly suggest that the possible values of the vortex numbers s± 
should depend on the spin-wave numbers r± by the following rule, 

• for r + + r_ odd — * s± = 0, ±1, ±2, . . . 

13 5 

• for r + + r_ even — ► s± = ±-, ±-, ±-, . . . (46) 

In order to investigate the validity of the proposal (I45|46l) beyond the decoupling point, we 
have solved numerically the original Bethe ansatz equations ( 1241) for L ~ 24. This numerical 
work enables us to compute the sequence, 

that are the expected to extrapolate to the dimensions Xy+'*z{y)- 

In table 1 we show the finite size sequences ( 1471) for some of the lowest dimensions with 
r + + r_ = 1 for 7 = |,|. The data for ^^'1(7) is restricted to L = 16 due to numerical 
instabilities with the respective Bethe roots. 
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L 


*o°'i°(f) 


<i(f) 




<i°(f) 






8 


0.272854 


1.449160 


1.494484(*) 


0.267351 


1.393699 


1.548692(*) 


10 


0.269116 


1.504508 


1.510610(*) 


0.263877 


1.434436 


1.568194(*) 


12 


0.266716 


1.542057 


1.518366(*) 


0.261697 


1.460574 


1.578323(*) 


14 


0.265022 


1.569286 


1.522338(*) 


0.260192 


1.478663 


1.584066(*) 


16 


0.263751 


1.590025 


1.524411(*) 


0.259085 


1.491841 


1.588751(*) 


18 


0.262754 


1.606393 





0.258233 


1.501900 





20 


0.261945 


1.619746 





0.257555 


1.509764 





22 


0.261272 


1.630840 




0.256999 


1.5161128 




24 


0.260702 


1.640206 




0.256535 


1.521357 




Extrap. 


0.250(1) 


1.811(2) 


1.52(1) 


0.250 0(1) 


1.581(1) 


1.59(1) 


Exact 


0.25 


1.812 5 


1.5 


0.25 


1.583 3... 


1.583 3... 



Table 1: Finite size sequences for the extrapolation of anomalous dimensions of the U q [Osp{2\2)\ model for 
7 = 7r/5, 7r/4. The exact expected conformal dimensions are -^0(7) = 1/4, ^Jifl) = 1/4+1/[4(7/tt)(1 — 7/77)] 
and -X"o'i(7) = 1/4+ 1/(1 — 7/77). The symbol (*) refers to Lanczos numerical data. 

In Figures 4(a, b) we exhibit the pattern of the roots associated to the X®'® (7) and X 'l (7) 
respectively. 
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In table 2 we show similar results for dimensions where r + +r_ is even and the corresponding 
roots structure are exhibited in Figures 5 (a, b, c). 



L 


*!,o 5 (|) 




^iV(f) 


*o7(f) 






8 


0.320662 


1.083548 


1.459096 


0.339496 


1.060728 


1.394641 


10 


0.318225 


1.092890 


1.530066 


0.337373 


1.068928 


1.452543 


12 


0.316831 


1.098146 


1.580667 


0.336207 


1.073385 


1.492308 


14 


0.315946 


1.101419 


1.618909 


0.335494 


1.076069 


1.521388 


16 


0.315341 


1.103607 


1.649054 


0.335025 


1.077806 


1.543652 


18 


0.314906 


1.105151 


1.673578 


0.334699 


1.078993 


1.561296 


20 


0.314580 


1.106286 


1.694025 


0.334462 


1.079839 


1.575661 


22 


0.314327 


1.107148 


1.711407 


0.334284 


1.080463 


1.587612 


24 


0.314126 


1.107821 


1.726419 


0.334147 


1.080936 


1.597730 


Extrap. 


0.312 53(1) 


1.112 3(2) 


2.054(1) 


0.333 34(2) 


1.833 2(2) 


1.752(1) 


Exact 


0.312 5 


1.112 5 


2.05 


0.333 3... 


1.083 3... 


1.75 



Table 2: Finite size sequences for the extrapolation of the anomalous dimensions of the U q [Osp(2\2)] model 

ii i i 

for 7 = 7r/5,7r/4. The exact expected conformal dimensions are X 2 q 2 (7) = 1/[4(1 — 7/71")], X^\ 2 (j) 

. 1 - 7 /tt) + 1/[4(1 - 7/tt)] and X*f * (7) = (1 - 1 /ix) + l/(4 7 /7r). 

All these numerical results confirm the conjecture f)45|46l) for the finite size properties of the 
U q [osp(2\2)] quantum spin chain. We now proceed with a discussion of our results. For periodic 
boundary conditions the ground state E^\l,j) sits in the sectors r + = ±1 and r_ = or 
r + = and r_ = ± and according to the rule (j4"6]) the respective vortex numbers have the 
lowest possible value s± = 0. From Eqs. (I44II45P we derive that its finite size behaviour is, 

= eg (7) + ^i(7) + O(L-). (48) 

By comparing Eqs. (137)480 we conclude that, in the continuum, the U q [osp(2\2)} vertex model 
should be described in terms of a conformal field theory with central charge c = —1. The 
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respective dimensions of the primary operators depends on the anisotropy and measuring 
them from the ground state we find that they are X£f'rz(l) ~ \ where r± and s± satisfy the 
condition (146j) . This is probably the first example in the literature of a theory with c < 
exhibiting a line of continuously varying exponents. In particular, we see the lowest conformal 
dimension occurs in the sector r + = r_ = with value X\ = ^- (1 Zi) which degenerates to 
that of the ground state for 7 = 0. The isotropic point 7 = possesses indeed special features. 
From Eq. (|45p we see that for s + 7^ s_ the scaling dimensions diverge as 7 —>■ 0. In this limit 
one expects therefore that only the sectors s + = s_ will contribute to the low-energy operator 
content. In order to describe the expected scaling dimensions at the isotropic point lets us, by 
considering the rule (j4"B|) . define s + + s_ = 2s (2s + 1) for r + + r_ = 2f + 1 (2f). We see then 
that the finite part of the dimensions (T4"5l) becomes, 



where f, s = 0, ±1, ±2, . . .. 

The above conclusions for 7 = agree with only part of the recent predictions made in 
[T4] for the possible values of the conformal dimensions of the isotropic osp(2\2) spin chain. 
Although the dimensions (1491) are the same for both sectors and coincide with that of a free 

n 

boson with radius of compactification R = 1 or R = 2 □, the respective values for the spin- wave 
or vortex numbers are restricted solely to odd integers. This subtlety may be of relevance in 
the description of the continuum limit of the osp(2\2) spin chain [31J. 



equations (T2"4"j) are not fully decoupled. We shall therefore start our study by considering the 
analytical predictions that can be made within the density root method. Considering that 
such an approach has already been well described before [28l [29] . we will present here only the 
2 Recall that the conformal dimensions of a compactified free boson (f>(x) — (f>(x) + 2-kR are " L -^ L - + ^ [SO]- 




2f + 1 



2f 



(49) 



5.2 The U q [osp(2\2m)} model 



For arbitrary m > 2 the special point 7 = f 



is not of great help because the Bethe ansatz 



19 



final results for general m. This framework can be adapted to handle the nested form of the 
Bethe equations (|2^|) and the finite size corrections for E^ m '(L, 7) turns out to be, 

(m+ 1) 



^)=«S>(7) + £*.<7> 



12 



+ X. 



Sl,...,S m _l,S + ,S_ 



l )•••)' m— ij* 



(7) 



+ L 



-2 



(50) 



where the corresponding scaling dimensions are given by 

-1 m— 1,± m— 1,± 



A ri,...,r m _ 1 ,r + ,r_ 



(7) = 7 E ^cS(7)^+ E *« ^ M (7) 



a,/3=l 



and the non-null matrix elements C^r (7) are 

7x r 7 



a,/3=l 



(51) 



07) 



2(1 - -)5 Q ,/3 - (1 - -) [<L,/3+l + ^,/3-lJ 
7T 7T 

7 

— (1 ) a = m — 1, /3 = ± or a = ±, /3 = m — 1 

* (52) 
1 a = / g = ± 



a, /3 



1, . . . ,m 



'1 



7T 



a 



±, /5 = T 



Taking into account previous experience with the m = 1 case, one would expect the exis- 



tence of some rule relating the possible values of the sets {si, . . . , s m -i, s+, s_} and {ri, 



-1,? 



For m > 2 we encounter some difficulties to unveil possible constraints between these numbers 
solely on basis of numerical solutions of the Bethe ansatz equations (|24|) and exact diagonal- 
ization of the respective Hamiltonian (!25l) . In order to make some progress we assume that 
the origin of such rule should go back to the issue of treating strictly periodic boundary con- 
ditions for the fermionic degrees of freedom in all sectors. We can first consider the situation 
in which all the 2(m + 1) degrees of freedom behave as bosons as far as boundary conditions 
are concerned. Next we look at the sectors whose eigenenergies do change as compared with 
the original system containing two bosonic and 2m fermionic degrees of freedom. The sectors 
whose spectrum remain the same should be described by integers while the remaining ones by 
half-integers as far as the values of s a are concerned. Having in mind the above considerations, 
we are able to derive the following conjecture for the constraints 



for r i+ i + rj_i odd 
for r i+ i + rj_i even 



Si 



1 3 5 
±-,±-,±-. 

2' 2' 2 

0,±1,±2, .. 



m — 2 
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• for r m _ 2 + r + + r_ odd 

• for r m _2 + r + + r_ even 

• for r m _i + r + + r_ odd 

• for r m _! + r + + r_ even 

Before proceeding we would like to note that the above constraints reflect the structure of 
the Bethe ansatz equations (I24p . In fact, the vortex numbers s a depend on the values of the 
neighboring spin-wave numbers according to the Dynkin diagram of Figure 3. Furthermore, 
such relationship for the Bethe roots with (O) o r without (®) self-scattering are just the 
opposite. In order to give some support to this conjecture we solve numerically the Bethe 
ansatz equations ( |24|) in the cases m = 2, 3 for some of the low- lying energies. In table 3 we 
have presented the results of the extrapolation for three possible dimensions for m — 2. In 
table 4 similar data is shown for m = 3. 



1 3 5 

S m _! - ±-,±-,±-, . 

s TO _i = 0,±l,±2, ... 

s± = 0,±l,±2, ... 

s+ = ±~, ±-, ±— , . . . 
± 2 2 2 



(53) 
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L 


A l,0,0l5 ) 


nil 

~Y U ' 2' 2 (ir\ 

^0,0,0 KJJ 


V 2 ' 2 ' 2 I TT \ 

A 1,0,1 I5J 


x-0,0,0/7r\ 
A l,0,0l4^ 


n — — 

X^"' 2 ' 2 (TT\ 
^0,0,0 \J) 


V 2 ' 2 ' 2 I 7T \ 

A i,o,i \T) 


8 


0.400274 


0.587137 


0.646331 


0.375903 


0.642602 


0.650216 


10 


0.400271 


0.594580 


0.635760 


0.375660 


0.648899 


0.640274 


12 


0.400230 


0.599575 


0.628254 


0.375490 


0.652331 


0.633314 


14 


0.400192 


0.603158 


0.622594 


0.375376 


0.655429 


0.628133 


16 


0.400160 


0.605863 


0.618139 


0.375296 


0.657300 


0.624103 


18 


0.400136 


0.607968 


0.614520 


0.375239 


0.658686 


0.620865 


20 


0.400116 


0.609660 


0.6115071 


0.375196 


0.659768 


0.618197 


22 


0.400100 


0.611046 


0.608950 


0.375165 


0.660590 


0.615954 


24 


0.400088 


0.612220 


0.606745 


0.375140 


0.661270 


0.614038 


Extrap. 


0.400 04(1) 


0.623(3) 


0.562 2(3) 


0.375 1(2) 


0.664(3) 


0.583 1(3) 


Exact 


0.4 


0.625 


0.562 5 


0.375 


0.666 6... 


0.583 3... 



Table 3: Finite size sequences for the extrapolation of anomalous dimensions of the ?7 g [Osp(2|4)] model for 



7 = 7r/5, 7r/4. The exact expected conformal dimensions are 

x i o(7) = (!-7/7r)/2, Xq.oo 2 (7) = 1/[2(1-7A)] 



1 1 



and Xjf'tfr) - 1/4 + 1/[4(1 - 7 /tt)]. 



All of them are in accordance with that predicted by Eqs. (I51II52P provided the rule (|53|) is 
taking into account. 

As before, we do not expect that the ground state Eq^^L,^) for general m will lie in the 
sector with all null spin-wave numbers due to the constraints ( 153)) . For m = 2, combination 
between Bethe ansatz and exact diagonalization results leads us to conclude that the ground 
state sits indeed in the sectors T\ = ±1 and r± — or t\ — r + — r_ = ±1. We have verified, 
for instance, that the lowest energy in sectors {1,0,0} and {1,1,1} are exact the same for 
finite L. Consequently, from Eqs. (1 5 |l 5 1 j) we derive that Eq (£,7) behaves as, 



£ {2) (£,7) 



3 (m) 



(7) 



7T 



■v 2 {i) 



3 - 6(1 - ^) 

7T 



(54) 



L ~°° " !/ GL 2 

We see that the term in Eq. (154"]) . usually related to the central charge, now varies continu- 
ously with the anisotropy 7. This is the typical expected behaviour for the critical properties of 
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L 


V 0,0,0,0/7T\ 

^•1,2,1,2 {J J 


n I L _ L 

V ' 2 ' 2 ' 2 / 7T \ 
A l,2,l,l {5) 


-1^0,0,0,0 /7T\ 

^1,2,0,3 I 5 ) 


-^0,0,0,0 /tt\ 
^1,2,1,2 I 4 ) 


— —— — — 

V u ' 2 ' 2 ' 2 / ir\ 

A l,2,l,l I4J 


-^0,0,0,0/7^ 
^1,2,0,3 I 4 ) 


8 


0.637377 


0.682230 


1.069574 


0.609176 


0.689721 


1.109450 


10 


0.643170 


0.700616 


1.067823 


0.615581 


0.697068 


1.115921 


12 


0.645843 


0.704425 


1.065523 


0.618706 


0.700729 


1.119000 


14 


0.647272 


0.706610 


1.063470 


0.620474 


0.702871 


1.120713 


16 


0.648119 


0.708095 


1.061765 


0.621580 


0.704113 


1.121774 


18 


0.648658 


0.708884 


1.060367 


0.622321 


0.705086 


1.122481 


20 


0.649018 


0.709568 


1.059214 


0.622843 


0.705583 


1.122975 


22 


0.649269 


0.710195 


1.058257 


0.623225 


0.706178 


1.123337 


24 


0.649448 


0.710345 


1.057453 


0.623513 


0.706353 


1.123608 


Extrap. 


0.650 05(1) 


0.711 8(3) 


1.052 (1) 


0.624 9(2) 


0.708 0(2) 


1.125 2(1) 


Exact 


0.65 


0.712 5 


1.05 


0.625 


0.708 3... 


1.125 



Table 4: Finite size sequences for the extrapolation of anomalous dimensions of the ^[05^(216)] model for 
7 = 7r/5,7r/4. The exact expected conformal dimensions are -X^'i'^ (7) = (3 — 2~f/ir)/4, 1 2 ' 2 (7) — 

(1 - 7 /tt)/2 + 1/[4(1 - 7/tt)] and = (3 + 6 7 /tt)/4. 

loop models [32] derived from vertex models with appropriate boundary conditions such as the 
g-state Potts and six- vertex systems [33]. The criticality of the loop model depends on its fu- 
gacity per every loop which turns out to be a function of the anisotropy 7 of the corresponding 
vertex model, see for examples [34J. In our case, strict periodic boundary conditions for both 
bosonic and fermionic degrees of freedom should probably work as the bridge between the loop 
and the vertex model formulations. This is at least the situation of the isotropic osp(2\2m) 
vertex model which was shown to provide a realization of an intersecting loop model with 
fugacity Q = 2(1 — m) [12J. Lets us admit that this analogy could be in some manner be 
extended for arbitrary 7. Considering that the U q [osp(2\2m)] vertex models share a common 
underlying braid-monoid algebra [17] it is natural expect that respective loop fugacity will be 
a function of the weight Q of the monoid operator. From our previous work [T7] it follows that 
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such weight is Q = —2 sin [7(771 — 1)] • Therefore, it is only at the special case m = 1 

that Q does not depend on the anisotropy, explaining why in this case the central charge was 
indeed independent of 7. A more precise description of these loop models such as the relation 
between vertex and loop Boltzmann weights has eluded us so far. 

We have carried out the above analysis up to m = 3. This leads us to conjecture that for 
general m the finite size correction for the ground state will be, 



m 1 

m+l -3(?n- 2[— ]- 

2 7T' 



(55) 



«7)-^n(7) 

where [y] denotes the largest integer less than y. 

We note that the result (|55|) when 7 — » agrees with the central charge behaviour predicted 
in [12] for the isotropic osp(n\2m) spin chains. In this limit, we also see from Eqs. (I51II52I) 
that the scaling dimensions for s + 7^ s_ diverge as 7 approaches zero and as before only the 
sectors s + = s_ contributes to the low- lying operator content. The generality of this scenario 
for arbitrary m strongly suggests that the continuum limit of the osp(2\2m) spin chains should 
be described by some peculiar field theory. In fact, we remark that a proposal towards this 
direction have recently been put forward in the work 



6 Conclusions 

In this paper we have studied an integrable vertex model invariant relative to the U q [osp(2\2m)] 
quantum superalgebra. The corresponding transfer matrix eigenvalue problem has been solved 
by the algebraic Bethe ansatz for a variety of grading choices. We thus have complemented 
previous efforts concerning the exact solution of solvable vertex models based on super algebras. 

We have explored the results for the transfer matrix eigenvalues and Bethe ansatz equations 
to investigate the thermodynamic limit properties as well as the finite size corrections to the 
spectrum in the massless regime. We have argued that the root density method needs a 
subtle adaptation to predict the correct finite size effects. It was observed that the constraints 
between spin-wave and vortex numbers are reflected in the Dynkin representation of the Bethe 
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ansatz equations. We believe that this will be the general scenario for integrable models based 
on superalgebras. This analysis has been helpful to point out possible classes of universality 
governing the criticality of the massless phase. The continuum limit of the £7 g [osp(2|2)] vertex 
model appears to be described by a c = — 1 conformal theory with critical exponents varying 
with the anisotropy. On the other hand, the gapless regime of the U q [osp(2\2m)] models 
for m > 2 was found to have a multicritical behaviour typical of loop models of statistical 
mechanics. 

We hope that our results will open further possibilities of investigations. For instance, one 
could use the Bethe ansatz equations to study the free-energy thermodynamics of the 
Ug[osp(2\2m)] vertex models. This representation is in fact rather suitable for the application 
of the so-called quantum transfer matrix method for finite temperatures [351 36J . This would 
provide us information on relevant physical properties of the U q [osp(2\2m)] spin chains such as 
specific heat and magnetic susceptibility in the entire temperature range. In particular, this 
could be used to check the S'-matrix of a osp{2\2) field theory proposed to described certain 
disordered systems [37J. 
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Appendix A: Two- parameter U q [osp(2\2)} vertex model 



In this appendix we present the algebraic Bethe ansatz solution of a two-parameter osp(2\2) 
vertex model. These parameters are directly related to the q- deformation and to the continuous 
U(l) parameter of the four-dimensional representation. The respective ^-matrix in the FBBF 
grading can be written as follows, 

= ^2a a {\)e aa (g) e aa + ^ d a ^(X)e^ a ^ 



4 4 

, p <8> e Q: 5 -p 

a=l a,P=l 

+ 61 (A) (I12 <S> e 2 i + e 2 i <8> e i2 + e 24 <8> e 42 + e 42 <g> e 24 ) 



+ b 2 (X) (ei 3 ® e 3 i + e 3 i <g> ei 3 + e 34 <S> e 43 + e 43 ® e 34 ) 

+ Ci(A) (e 22 <g> en + e 44 <g> e 22 ) + c 2 (A) (e 33 <g> e n + e 44 <g> e 33 ) 

+ ci(A) (en ® e 22 + e 22 ® e 44 ) + c 2 (A) (e n (g) e 33 + e 33 ® e 44 ) 



(A.I) 



The main Boltzmann weights are given by 
0l (A) = -i (e 2A g 2 - l) (e^-l) 
03(A) = 1 (e 2 \ 2 - l) (e 2A - g 2 ) 
W = I (e 2 * - l) (e 2A g 2 -l) 
ci(A) = -\ (e 2 V - l) (g 2 -l) 



2A, 



Ci(A) = e 2A Cl (A) 
while the remaining elements can be written as 



02(A) : 

a 4 (A) 

62(A) 

c 2 (A) 
c 2 (A) 



yj) 



e 2 \l - l) 



^( e2 V"l)( 

l(e--l)(e 2 ^ 2 -l) 
-1 (e 2A g 2 - l) (ql - l) 



e 2A c 2 (A) 



(A.2) 



-2i(e 2A 

91 V 



£ (e 2A - l) (e 2 ^ 2 - l) 
-i{^-l)(l-q^f{l-q^ 
g 2 e 2A (e 2A -l) (l - q?) h (1 - q^ 

-^(g 2 -i)(? 2 2 -i) 



1 



1 **(«?- 1) («2 
e 2A (1 - qjql) + qf + q\ - 2 

e 2A {ql + ql - 2q\q 2 2 ) + o 2 g 2 - 1 



a = (3 = 1,4 
a = /3 = 2,3 

a</3, /3 - a = 1,2 
a > /3, a - /3 = 1,2 
= 5-/3 = 2,3 
a = 5 - /3 = 1 
a = 5-/3 = 4 



(A.3) 
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It is not difficult to see that for qi = g 2 = Q one recovers the U q [osp{2\2)} R- matrix defined 
in Eqs. (|3]17j) when the grading FBBF is adopted. We start by recalling the definition of the 
monodromy operator entering the algebraic Bethe ansatz solution of this vertex model in the 
presence of inhomogeneities, 



TO(A, fa}) = r { {1(\ - ^)i?!£_!(A - /^_0 . . . R { {!(\ - a*x), 



(/) 



(A.4) 



as well as the associate row-to-row transfer matrix 



rO(A,fo}) = Str [TO (A, {/,,}) 



(A.5) 



The monodromy operator ()A.4j) plays an important role in the formulation of the algebraic 
Bethe ansatz method, and with the help of the Yang-Baxter equation one can show that it 
satisfies the following quadratic algebra 



R{{\\ - /x)T^(A, ® T (/) (//, {„,-}) = TW( ft {/*,}) £. T^(A, {^})R[ f 2 \\ - n). (A.6) 



CO, 



We remark that the super tensor products in (1A.6j) takes into account the parities in the 
grading FBBF [19]. Besides that, another important ingredient for an algebraic Bethe ansatz 
solution, is the existence of a pseudovacuum state |$o) i n which the monodromy matrix acts 
triangularly. For the considered vertex model (IA.lHA.3l) we can choose 

/1\ 



l < fo)=(g)|0) J 



in which the action of the operator Rj^(X) gives 



fuiM |o) j 







t 

w»(A) |0>, 











t 



«*W 10), 





t 
t 

t 

c*(A)|0)J 



(A.7) 



(A.i 
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The symbol f stands for non-null values while the functions uJ a {X) are given by 

^i(A) = - ai (A) w 2 (A) = 61(A) 

u 3 (X) = 63(A) o; 4 (A) = -d AA (X) (A.9) 

Previous experience with similar vertex models [16] leads us to adopt the following repre- 
sentation for the monodromy matrix flA.4j) 

( B(X, fa}) B 1 (X, fa}) B 2 (X, fa}) F(X, fa}) \ 
C^fa}) A n (\fa}) A 12 (X,fa}) B*(X,fa}) 
C 2 (X,fa}) A 21 (\fa}) A 22 (\fa}) B*(\fa}) 
V C(X, fa}) C?(A, fa}) C*(A, {^}) D(X, fa}) J 
and the diagonalization problem for the transfer matrix becomes equivalent to the problem, 



T {f \X,fa}) 



(A.10) 



-B(X, fa}) + E At (A, W) - £>(A, {^}) 



8=1 



A (/) (A,{^}) 



(All) 



The triangular form exhibited by flA. 8|) together with (IA.4I) allow us to compute the action 
of elements of the monodromy matrix T^(A, {/%}) on the pseudovaccum state |$o)- In this 
way we can regard Bi(X, {f-ij}), B 2 (X, {fij}) and F(X, {fJ-j}) as creation fields while the diagonal 
ones satisfy the relations 



B{\,fa})\*o) 



L L 

IJ wi(A - //,) |$ 0) -D(A, {^}) |$ ) = II ^(A - fa) |$ ) 



1=1 

L 



i=l 



Ai(A,{//j}) |$o) = n>2(A-//i) |$o) A 22 (A,{/i i })|$ Q ) = n^3(A-^)|$ ), 



8=1 



«=i 



(A.12) 



as well as annihilation properties for the remaining elements 

C(\,fa})\$ ) = A lJ (X,fa})\%) = i^j 

C t (X,fa})\$ ) = C*(X,fa})\® ) = (A.13) 

The above relations imply that |$o) is an eigenstate of the transfer matrix whose respective 
eigenvalue is 

L L L L 

A? (A) = - I] u;i(A - m) + II ^(X -&) + ![ u 3 (X -fn)-U cu 4 {X - m). (A.14) 



i=i 



i=i 



i=l 



i=l 
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Within the algebraic Bethe ansatz method we now look for the remaining transfer matrix 
eigenvectors as linear combinations of products of creation fields acting on |$o)- The general 
form of these eigenvectors has been already presented in [16]. In order to accomplish that we 
need to disentangle from Yang-Baxter algebra flA.61) appropriate commutation rules between 
the diagonal and creation fields. Until this stage, this approach is quite similar to the one used 
in [16] unless by the fact that the presence of two deformation parameters modifies the set 
of commutation rules required. In order to avoid an overcrowded section, these commutation 
rules have been collected in appendix B. 

A careful analysis of the commutation rules given in appendix B, together with the relations 
(IA.12IIA.13IIA.14I) . leave us with 

A(/)(A,{^}) 



- [] wx(A - in) II - A ) a b(A, {A,}) + A A {\, {A,}) 

1 



t=i 

L 



i=l 



-J]o; 4 (A-^)n 



A D (A,{A,}) 



(A.15) 



»=1 i=l ^4,4 (A - A 4 ) 

where A#(A, {Aj}), A D (A, {Aj}) and A^(A, {Aj}) are eigenvalues of the auxiliary matrices T B (X, {Aj}), 
To (A, {Aj}) and Ta(A, {Aj}) respectively. The set of rapidities {A.,} follows from the vanishing 
condition of the so called unwanted terms which will be discussed later. 

Initially we shall consider the auxiliary transfer matrix T^(A, {Aj}) defined as 



T A {\ {\}) = Tr [G(X, {A l })r Ql (A - A 1 )r a2 (A - A 2 ) . . . r an (A - A n )] 
with the following structure for the auxiliary r-matrix 



r(A) 



/a? (A) 











b* (A) 











6*(A) 


V o 








rtits are 







o \ 




2(A)/ 



a{(\) 



yi 



(A.16) 



(A.17) 
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4(A) = 
and the diagonal twist is given by 



0( 
1)( 



e 2A - ft 2 



e 2 \ 2 2 - 1 



(A.18) 



G(A,{A,}) 



V 



nf=i^ 2 (A-^) 

nr=i &i(a-a 4 ) 







nf=i ^3 (a - 



(A.19) 



nti& 2 (A-A, 

We are interested in the solution of the eigenvalue problem 

T A (X,{X t })f = A A (X,{X t })f 



(A.20) 



which is trivial due to the diagonal form of flA.17[) . 

Defining the spin up state |T) — I I an d the spin down state ||) = I | wo can write 



•F=iT)®IT>®---U>®---®IT> 



(A.21) 



possessing n+ spin up states and n_ spin down states such that n = n + + n_. In this basis 
it is convenient to separate the set of rapidities {Xj} into two subsets {A^} and {Xj}, each 
one associated with the spin up and spin down components of T respectively. With the above 
considerations we have 



TT (\ \ TT ^*(A ~ Aj ) -r-r ^(A - Aj ) 



(A.22) 



l 1 = 1 1 &2(A-A/-) l t L 1 6 2 (A-Ar) 

Next we turn to the auxiliary matrices Tg(A, {Aj}) and Tb(A, {Aj}), and their respective 
eigenvalues. By way of contrast, these matrices are not defined as a trace of a monodromy 
matrix, but they are diagonal matrices whose elements are given by 

1 



T B {\{\mizz = n 



...5, 



OtnPn 



n 

T D (X, {Aj})^ 2 '''!™ = J| b ai (A - Xi)5 ail3l 5 a2 f3 2 . . . 6 anf3n 



(A.23) 



i=l 
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Considering the trivial eigenvectors T (1A.21I) we are left with 

1^1 



A D (A,{A,}) 



n 



n 



n+ n_ 



(A.24) 



i=l i=l 

In this algebraic Bethe ansatz construction the unwanted terms are canceled out by making 
use of explicit form for T (IA.2ip and provided that the set of rapidities {Xf} and {A^} satisfy 
suitable Bethe ansatz equations. Putting our results together, the eigenvalues A^(A, {fa}) are 



A (/) (A,M) 



- J]o;i(A-/i 

L 

i=i 

L 

+ W^X-fa) 

i=l 

L 

- n w 4(A-/ij) 



A-z'f g_ (A -if 




q_ (a + z 7 2 + if) 



X + if) Q_ A + if 



A + z 7l +zf) g_ (A -if) 



Q_ (A + 



71 



1=1 



(A + ryi + i 


72^ 

2 ; 


Q- 


(A + z'72 + i 


71 A 
2 y 


Q+ 


(A + zf) 


<2- 


[x + if) 





(A.25) 



n± 

where as in the main text Q±(A) = [] sinh (A - A| 



(±) 



i=l 



The corresponding Bethe ansatz equations for the variable {A^- } are given by 



n 



sinh fA+ - fa + if 



i sinh ( A 



n 

i=i 



sinh (A 



fa 



1 2 



n 



sinh (Xj 



A" + z 



(71+72) 
2 



, 1 sinh ( AJ - A, 



fa + if 



sinh [\j 



fa i\ 



n + sinh ( X • 

n 1 



A+ + z 



; (71+72) 
2 

(71+72) 



i=l 



sinh \ Xj 



Xt 



,■(71+72) 



(A.26) 



We finally remark that in the above relations we have set qi$ = e* 71,2 and considered the 
shifts {A+} — > {A+} — if and {Xj} — > {Xj} — if. In order to obtain the nested Bethe ansatz 
solution of the U q [osp(2\2m)] vertex models one has to use Eqs. flA.25IIA.26p in the final step of 
the recurrence relations flll|ll2l) at the particular point 71 = 72 = 7- 
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Appendix B 

In this appendix we have collected the set of commutation rules required to perform the 
algebraic Bethe ansatz for the two parameters U q [osp{2\2)] presented in appendix A. 



B {m{ ») = _g^^)B(A) + ^5,(A)B0x) (B.l) 



+ iB^7) mcKX) ~ %^^< A >^> (R2) 



+ 



(i 4 ,4(A-^) J ' * rf 4)4 (A-/i) 6j(A-/i) 

1 



6j(A - //) 



d 4j -+i(A - ^)d M (A - ^) 

" ^ MM 



^■ 3 - fe F(^a(A) 

(B.3) 



B^)F(A) 



D(A)ityx) 



+ 



rf 4 ,4(A - /i) 

- /i) 

<i 4j4 (A — /i) 
a 4 (A — /i) 



4,4 (A - /i) 
4,i (A - /i) 



F(^)S(A) 



c? 4)4 (A — /i) 

4,i+l(A - /i) 

4,4 (A - //) 



F(jm)D(X) 



4,4 (A 



F(\)D(p) 



(B.4) 



(B.5) 



^•(A)F(/i) = 



6-,(A — /x) Cj(X — //) Cj(A — /i) 



MA-Ai) 
q(A - /i) 



bi(X - /u) 6j(A - fi)_ 
Cj(X-n) 



bi(X - fj) 



b 3 {x)b;^)- 



bj(X-n) 



FiriMX) + 
B*(X)Bj(ri 



Cj(X - /i) cj(X - fj) 
k(X - n) bj(X - n) 



F{X)A^) 
(B.6) 
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a\(X — /i) 
1 

ai(A — /i) 



(A - /i) 



rf4,4(A — /i) 

^4,i+i(A - /x)rfi [4 (A - /^) 

^4,4(A — /i) 



S^F^BiX) 



(B.7) 



[F(A),F(/i)] = (B.8) 

F^)^(A) = - ^:^ t (A)^) + ^^S,(/i)F(A) (B.9) 

Bi(n)F(x) = - ^"^ (W) + ^x^y F ^^( A ) (B.10) 

+ ~ ^ (A)C7i(M) ~ " fa s-^(A04*(A) (B.11) 

04,4 (A — /ij "4,4 (A — /i) 

~ ( R12 > 
C ''W B » - -H3 Bj0 " a(A) - M3 B( " )Ai(A) 

C'tAJfl^) = 1 B t (/i)<?(A)RS(A - p) - f^Z^ BiWCm 

- ^?*« B ^(A) (B.14) 

In order to clarify our notation, the elements are obtained from flA. 17|) through the 
definition 

KA)= E r2(A)^®e«. (B.15) 

i,j,k,l 
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Finally, we have also used the relation TZ^X) = ^Xy+iM where the elements Rj£ t follows 
from the convention 

4 

R{\) = J2 Rfi{\)ei3®eu. (B.16) 

i,j,k,l 
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Figure 4: The Bethe ansatz roots Xy' (o) and A] ' (x) for 7 = | and L = 12. The roots refer to the 
dimensions (a) JTg'°(7) and (b) X^l(-j). 
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Figure 5: The Bethe ansatz roots ' (o) and \\ ' (x) for 7 = | and L = 12. The roots refer to the 
dimensions (a) (7) , (b) X? (7) and X?y* . 
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